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We calculate the interaction potential between two colloids immersed in an aqueous mixture 
containing salt near or above the critical temperature. We find an attractive interaction far from 
the coexistence curve due to the combination of preferential solvent adsorption at the colloids' 
surface and preferential ion solvation. We show that the ion-specific interaction strongly depends 
on the amount of salt added as well as on the mixture composition. Our results are in accord with 
recent experiments. For a highly antagonistic salt of hydrophilic anions and hydrophobic cations, 
a repulsive interaction at an intermediate inter-colloid distance is predicted even though both the 
electrostatic and adsorption forces alone are attractive. 



I. INTRODUCTION 

Charged surfaces in liquid media are ubiquitous in soft 
matter and the interaction between such surfaces has 
been studied extensively. In the seminal theory of Der- 
jaguin, Landau, Verwey, and Overbeek (DLVO) the total 
interaction potential is the additive combination of the 
attractive van der Waals interaction and an electrostatic 
repulsive screened Coulomb potential [HE]. In the DLVO 
theory the liquid between the surfaces is a homogeneous 
dielectric medium and only enters the electrostatic po- 
tential through its permittivity s. The situation is differ- 
ent in solvent mixtures, where the structure of the solvent 
may vary in space due to gradients in electric field and ion 
density. Several authors explored experimentally the in- 
teraction between charged surfaces in mixtures [5HZ]. In 
a homogeneous binary mixture of water and 2,6-lutidine 
close to the demixing temperature, a reversible floccu- 
lation of suspended colloids has been observed [3]. The 
flocculated region in the phase diagram is suppressed by 
the addition of salt [4]. These experiments were inter- 
preted in terms of the preferential adsorption of one of 
the solvents on the colloid surface in accord with wetting 
theory [8l |9] . More recently, it was suggested that the se- 
lective solvation of ions in the two solvents is important 
in these experiments [TUl E] as well as charge-regulation 
effects Q3] • 

Recently, direct measurement of the interaction poten- 
tial between a colloid and a charged wall was performed 
in water-2,6-lutidine mixtures below the lower critical 
solution temperature (LCST) [7]. In these experi- 
ments the addition of salt greatly influenced the range 
and strength of the interaction [TJ, leading to an attrac- 
tive interaction far from the critical temperature T c . Ad- 
dition of salt has a remarkable effect on the flocculation of 
colloids in a critical mixture of water-3-methylpyridinc, 
where it reduces dramatically the onset temperature of 
the flocculation [S] . Preliminary experiments |12j suggest 
that the variation of the mixture composition also has a 
pronounced effect on the interaction. 

A theoretical effort has been made to incorporate ionic 
solvation effects in polar mixtures [lOj [TTJ [13HT5] . The 
solvation energy of an ion arising from the ion-dipole in- 



teraction and sometimes other specific interactions with 
the polar solvent strongly depends on the solvent and 
the chemical nature of the ion [TBHT5] . The total solva- 
tion energy of an ion when it is transfered from one sol- 
vent into another (the Gibbs transfer energy) is in many 
cases much larger than the thermal energy, especially if 
one solvent is water and the other is an organic solvent 
[T?lU9| . In addition, for a given solvent, the solvation 
energy of cations and anions is generally different due to 
their different size and specific chemical interactions with 
the solvent, such as Lewis acid-base interaction. Thus, 
the magnitude and sometimes even the sign of the Gibbs 
transfer energy for cations and anions may differ in var- 
ious salts and solvents [TH1 HO] • 

In bulk binary mixtures, the addition of salt modifies 
the coexistence curve [2"TH2"5] . When a salt-containing 
mixture is confined between charged surfaces, the influ- 
ence of ionic solvation on the surface interaction in mis- 
cible solvents has been studied in Ref. Q21 while for par- 
tially miscible mixtures close to the coexistence curve it 
was explored in Ref. 1101 In the latter case it was shown 
that preferential solvation alone can qualitatively modify 
the interaction between similarly charged surfaces, re- 
sulting in a strong long-range attraction. 

In order to account for the recent experimental findings 
in critical mixtures [TJ, Bier et. al. [13] had employed a 
linear theory that includes both ion solvation and solvent 
adsorption on the confining surfaces. Recently, Okamoto 
and Onuki [TT] also included in the theory solvent de- 
pendent charge-regulation, and studied the interaction 
induced by prewetting and wetting transitions close to 
the coexistence curve. 

In this paper we theoretically investigate how the inter- 
play between ion solvation and solvent adsorption affects 
the interaction between surfaces in the nonlinear regime. 
We show that a strong attraction exists even far from 
T c due to the nonlinear coupling between solvation and 
adsorption. These effects are strongly ion-specific and 
hence depend on the salt properties and concentration. 
We highlight the qualitative features of the theory and 
give several new predictions. 

We focus on the regime far from the prewetting and 
wetting transitions and neglect charge-regulation, assum- 
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ing large surface ionization. Thus, we do not expect any 
discontinuities in physical quantities such as pressure and 
adsorption. We concentrate on electrostatics and ion- 
solvation forces and hence ignore the van der Waals in- 
teraction (which is relatively small in the current range 
of parameters). 

The organization of this paper is as follows. In Sec. |TT] 
we present a coarse-grained model of a salt-containing bi- 
nary mixture, taking into account the specific ion-solvent 
and surface-solvent interactions. In Sec. |III| we present 
results for critical and off-critical mixture compositions 
and show how ion specific effects can be used to tune the 
inter-colloidal potential. Conclusions are given in Sec. 
W 



II. MODEL 

The surfaces of the two colloids are modeled as parallel 
plates of area S separated by a distance D and carrying a 
uniform charge density ea^ and ecru per unit area, where 
e is the elementary charge and the indices L and R denote 
the left and right plates, respectively. An aqueous binary 
mixture containing salt is confined between the plates. <f> 
is the volume fraction of the water (0 < <fi < 1), while 
the densities of the positive and negative ion are denoted 
by n + and n~ , respectively. Given their low density and 
small size, the volume fractions of the ions are neglected. 

For the fluid free energy we have Ft, = J ft,dr, where 
ft, is the bulk free energy density: [23J US] 
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Here the Boltzmann constant is set to unity and T is the 
thermal energy, vq = a 3 is the molecular volume of the 
solvent molecules where a is the linear dimension of the 
molecules. The first term in Eq. ([!]) is the mixture free 
energy [27]: 

v f m = 01og(0) + (1 - (^)log(l - <f>) + x<Kl ~ <t>), (2) 

where x ~ V^ 1 i s the Flory interaction parameter. In 
this symmetric model the critical composition is given 
by (j) c = 1/2. The square-gradient term in Eq. takes 
into account the free energy increase due to composition 
inhomogeneities [27]. where C is a positive constant. The 
third term in Eq. (n]) is the electrostatic energy density, 
where ip is the electrostatic potential. This energy de- 
pends on the composition through the constitutive rela- 
tion for the dielectric constant: e — e{<j>). For simplicity 
we use the linear form: s((f>) = e c + (e w — £ c )4>, where 
e w and e c are the water and cosolvent permittivities, re- 
spectively. The second line of Eq. ([I]) is the ideal-gas 
entropy of the ions, valid at low densities, while the third 
line is the ion solvation in the mixture. In the simple bi- 
linear form we employ, the strength of selective solvation 



is proportional to the coefficients Au + and Ait~ for posi- 
tive and negative ions, respectively |26j . In this form the 
Gibbs transfer energies of the ions are 
where </>i and </>2 are the compositions of two solvents. 

In addition to the bulk term, one adds the contribution 
of the confining surfaces to the free energy: F s — J f s dS, 
where f s is the surface free energy density given by: 



f s = eat/j s + Ajcj) s , 



(3) 



where ip s and <p s are the electrostatic potential and mix- 
ture composition at the surface, respectively. Eq. Q in- 
cludes a short range interaction between the surface and 
liquid which is linear in <£>, neglecting higher order terms. 
The surface field A7 measures the difference between the 
surface-water and surface-cosolvent surface tensions. 

The fluid is in contact with a matter reservoir having 
composition </>o and ion densities n . Hence, the appro- 
priate thermodynamic potential to extremize is the grand 
potential f2: 



n = F h + f s 



[fi<p + X + n + + A n ]dr, (4) 



where fi and A are the Lagrange multipliers for the com- 
position and ions, respectively. They can be identified as 
the chemical potentials of the species in the reservoir: 

i" = MoOfoj^o) and A± = Aq^O)^)' 

The governing equations are obtained by minimization 

of O, leading to the Euler-Lagrange equations 



+ logfuora* ) - Am 4> - A ± = 0, (5) 
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Au n — Au n — /x = 0. 



(7) 



The first equation yields the Boltzmann distribution for 
the ion density 



(8) 



Note the dependence on the mixture composition in this 
equation. The ion distributions are inserted into Eq. ^ 
to obtain a modified Poisson-Boltzmann equation with 
the boundary conditions — n • V?/»l,_r = gvl.r/z- n is 
the outward unit vector perpendicular to the surface. 
The governing equation for the composition [Eq. ([7|] 
is supplemented by the boundary conditions n- V<Pl,r = 
^1l,r- Assuming that ip — in the reservoir, the La- 
grange multipliers for the different species are: 



Mo 



df„ 



Au n, 







Au n 







Aq = logOonJ) - Ai 



(9) 
(10) 



In our effectively one dimensional system, solution of 
the Euler-Lagrange equations yields the density profile 
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4>{z) and potential ifj(z) from which we calculate all rele- 
vant quantities. Of particular interest is the pressure P n 
exerted on the surfaces by the liquid. — P n is the normal 
component of the Maxwell stress tensor and it is given 
by 01 |2g]: 



n (1 — Am 



f m +n+(l-Au+4>) 
(11) 



In a planar geometry P n is the zz component of the stress 
tensor and is uniform in mechanical equilibrium. The net 
pressure on the plates is given by the osmotic pressure 
II = P zz — Pb, where Pb is the bulk pressure. The inter- 
action potential between the plates U(D) is calculated 
from the osmotic pressure through 



U(D) =S Il(D')dD' 

J oo 



(12) 



The excess surface adsorption in a planar geometry is 
defined as usual by 



r = -^f [<p( z ) - <p ]dz. 



(13) 



This quantity is measurable experimentally and accounts 
for the total excess of fluid between the plates relative to 
the bulk. 



III. RESULTS 
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In the following we examine the interaction potential 
U(D) for a binary mixture containing lOmM of salt at 
temperatures larger than the critical temperature, T > 
T c . 

A. Colloidal interaction in a mixture at a critical 
composition 

We start with a reservoir at a critical composition, 
4>o = 4>c confined by two hydrophilic (Aj^.l > 0) col- 
loids. We first consider the symmetric case where both 
colloids are identical: A7L = A'f n and ox = or- We 
set a L ,R = -cr sa t, where cr sat = ^8n /(7rZ B ) and l B = 
e 2 /(Ane(^>o)T) is the Bjerrum length. er sat is used be- 
cause the effective surface charge saturates to this value 
in the high bare charge limit |29j . Furthermore, we as- 
sume symmetric preferential solvation, i.e. Au + = Au~ , 
and in this case the solvation asymmetry parameter, de- 
fined as Au d = Au + — Au~, vanishes: Au d = 0. 

The parameters we use are such that in the resulting 
profiles ip(z; D) and cf>(z; D) the conditions etp(z; D) /T > 
1 and Au (<Mz;D) — <j>o) > 1 hold. Hence, a lineariza- 
tion of Eq. ([5]) is not expected to be accurate [TTJ [T3] . 
In fact, we show below that it is the nonlinear coupling 



FIG. 1. (a) The interaction potential U(D) between two col- 
loids at a distance D at different temperatures r = T/T c — 1 > 
immersed in a mixture at a critical composition (<f>o — <j) c ). 
U(D) becomes attractive as r decreases. Here no = lOmM 
and A7b_£ = 0.1T/a 2 , corresponding to about 3.4 mN/m. 
For the solid curves, Ai^ = 4 and the surfaces have the 
same charge ctl,r = — o- sat . Dash-dot curve: the same as the 
solid curve for r = 0.008 except that <jl = 3o\r = — 1.5ow- 
Dashed curve: the same as for r = 0.008 except that Au~ = 
8. (b) The corresponding excess surface adsorption F. In 
this and in other figures, as an approximation of a water-2,6- 
lutidine mixture we used T c = 307.2K, v Q = 3.9 x 10" 29 m 3 , 
C = x/a [27], e 2 ,6-iutidino = 6.9 and e W ator = 79.5, and the 
surface area is taken to be S = O.Olum 2 . 



between the governing equations that is responsible to an 
attractive interaction far from T c . 

The solid curves of Fig. [lja) show U (D) as a function 
of the reduced temperature r = T/T c — 1 for T > T c . Far 
above T c U(D) is purely repulsive; as the temperature 
is decreased toward T c the interaction becomes attrac- 
tive. The change in the shape of U(D) originates from 
the interplay between surface field, the solvation induced 
attraction, and the electrostatic repulsion between the 
plates. The attraction between the plates is stronger 
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down to t — 0.003 and then becomes less attractive for 
t = 0.002. The attraction range of a few tenths of nm, 
its strength of a few T, as well as the temperature depen- 
dence are similar to those observed in recent experiments 

mm- 

The dashed and dash-dot curves in Fig. [l|a) are the 
same as the solid curve for r = 0.008 except for one pa- 
rameter. The effect of charge asymmetry is shown by 
the dash-dot curve, where the total charge was kept con- 
stant, but ctl = 3(7fj. Here the attraction is stronger, as 
in the classic Poisson-Boltzmann theory [30 , since the 
electrostatic repulsion between surfaces is weaker for the 
asymmetric case. The dashed curve shows that increas- 
ing Au~ by 4 (Au d = —4) results in a slightly weaker 
interaction. In the linear theory [TTJ [T3] , solvation cou- 
pling only enters as a term cx Au d . Thus, by comparison 
to the solid curve in Fig. [jja) where Au d = 0, it is clear 
that in the nonlinear regime this type of contribution is 
smaller than contributions proportional to Ait . 

The Gibbs adsorption V corresponding to the curves in 
Fig. [lja) is shown in Fig. [ljb). When the surface sepa- 
ration decreases T increases as the adsorbed fluid layers 
near the walls merge. Furthermore, since the density of 
the ions also increases, more fluid is drawn to the walls 
because of the solvation interaction. T increases as the 
adsorbed fluid layer thickness, which is comparable to 
the bulk correlation length, increases closer to T c . 

The parameters we use are close to those used in ex- 
periments in a salt containing water-2,6-lutidine mix- 
ture below the LCST [12]. In these experiments it 
is reasonable to assume that the surfaces have different 
charge densities because of their different chemical na- 
ture. Therefore, below we use ctl = 3<tr = — l-5c SQ t 
for hydrophilic surfaces. The value of A7 was chosen 
arbitrarily since there is no molecular theory to predict 
it accurately. Furthermore, in water-2,6-lutidine mix- 
tures the anions are expected to favor the water environ- 
ment more than the cations, since 2,6-lutidine is a Lewis 
base. This is supported by data of Gibbs transfer ener- 
gies of ions in water-pyridine mixtures [T5]. 2,6-lutidine 
is a structural analog of pyridine and hence on the basis 
of Ref. [Tniwe took Au + = 4 and Au~ = 8. We stress 
that over a wide range of values for <jl,r and positive 
Ajl,r, the results do not change qualitatively. This is 
also true for the solvation parameters as long as the ions 
are hydrophilic, Au* > 0. Note that the molecular vol- 
umes of the two components differ substantially in the 
real mixture leading to <j) c 7^ 1/2. Furthermore, our sim- 
plified model mixture has an upper critical solution tem- 
perature, whereas experiments are performed below the 
LCST in water-2,6-lutidine mixtures. Thus, one should 
interpret results in terms of the absolute distance from 
T c . 

U(D) can be attractive when only surface fields are 
present (Au* = 0, Ajl^ ^ 0) due to critical adsorp- 
tion [3lTf3"5] . When only preferential solvation exists 
(Au 7^ 0, A^l r — 0), we have shown recently [TUj that 
fluid enrichment near the surface, induced by solvation, 
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FIG. 2. The effect of removing preferential solvation or short- 
range chemical interactions on the potential U (D) between 
the colloids at temperatures given by (a) r = 0.009 and (b) 
r = 0.004. The dashed and dash-dot curves show U (D) when 
either the surface chemical affinity or preferential solvation 
parameters are zero, respectively. In the solid curves both 
short-range chemical preference and solvation are included. 
These two interactions are clearly non additive as the solid 
curve is not the sum of the dashed and dash-dotted lines. 



leads to an attractive interaction near T c . Interestingly, 
as is shown in Fig. [2j the attraction between the surfaces 
is significantly altered by the coupling of solvation and 
surface fields. 

In part (a) we highlight this idea by choosing param- 
eters such that far from T c , removing either the surface 
tension (dashed curve) or preferential solvation (dash- 
dot curve) the interaction is purely repulsive. When 
both forces are present the result is surprisingly a strong 
attractive interaction between the colloids (solid curve). 
Fig. [2^b) shows that closer to T c , an attractive interac- 
tion of similar magnitude exists when either force is miss- 
ing, but the interaction is unexpectedly much stronger 
when both are present. Thus, in the nonlinear regime, 
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FIG. 3. Interaction potentials at different temperatures r, 
(a) for <p = 0.48 < <j> c with A7 fliL = 0.1T/a 2 and (b) for 
cj>o = 0.52 > (j> c with A7fl,L = 0.4T/a 2 . The onset temper- 
ature for attraction is higher for <f>o < cj> c due to preferential 
solvation. In (a), at intermediate temperatures the potential 
has metastable states. Here we used Au + = 4, Au~ = 8 and 
o\l = 3ctr = — 1.5a S at- We took the average ion density to be 
no = lOmM leading to k, values of k ~ 2.69nm. 



solvation and adsorption forces are non additive and their 
coupling can account for attractive interactions far from 



B. Surface interaction in an off-critical mixture 

The interaction potential between two surfaces for off- 
critical mixture compositions at different temperatures 
is shown in Fig. [3| The legend also shows the di- 
mensionless parameter K^b corresponding to r, where 
K = ^/SttIbHq is the Debye wave number and is the 
bulk correlation length in the absence of ions, defined as 
6>(t) = ^fxa 2 / (d 2 f m (<l>c)/d<j) 2 ). The combination 
is useful for comparison with experimental results. Ex- 
perimentally, it was found that ions do not modify the 
correlation length significantly [7 . 
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FIG. 4. Solid curves show the dependence of U(D) on the 
mixture salt concentration no with Au + = 4 and Am" = 8; 
colloidal attraction increases with the addition of salt. The 
attraction is weak for a salt concentration of no = 0.01M but 
with Ati* = (dashed curve). Here cjf>o = (f> c , t = 0.008, 
A.jl,r = 0.1 and (Tl — 3<tr = —1.5a sa t- 



As is seen in Fig. [3] for O 0c the curve U(D) 
becomes attractive closer to T c , similar to the 4>o = 4> c 
case, but the temperature range at which this occurs is 
different. The attractive interaction is stronger for <po < 
4> c [Fig. [3] (a)] and thus the temperature range above T c 
at which it occurs is larger; compare with r values in Fig. 
[2] Note that in Fig. [3] (a) several curves have metastable 
states and that the interaction is repulsive at long range. 
The attractive interaction is weaker for 0o > 4>c [Fig. [3] 
(b)]. Here, we had to use a larger value of water-surface 
interaction in order to obtain a potential depth similar to 
the one which can be resolved in experiments (few T's) 

In order to understand the dependence on reservoir 
composition we insert Eq. ( 10 ) for Aq into the Boltzmann 
distribution Eq. ([8| to obtain 



n e 



=Fei/'/T+A u ± (0-0o) 



(14) 



This equation suggests that the ion density near the 
plates increases when 0o is reduced, which is verified nu- 
merically. This leads to a decrease of the osmotic pres- 
sure in the system since when Au ± > 1 the solvation 
contribution to n overcomes the ions entropic repulsion 
[cf. Eq. (11)]. Thus, the coupling between the mixture 
composition and the ion density is important for the in- 
teraction between the surfaces. Due to the ion-solvent 
coupling the attraction is increased when salt is added. 
This effect is shown by the solid curves in Fig. [4] and is 
seen in experiment [SJ Q2] . The dashed curve in Fig. [4] 
shows U(D) for the highest salt concentration but in the 
absence of ion-solvent coupling (Ait* =0). In this case 
the attraction is very weak, indicating that the decrease 
in Debye length when salt is added is not significant for 
the potential at the given parameters. 
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In general, attraction occurs up to a temperature win- 
dow AT = T - T c « 5K above T c , quite far from T c but 
not as large as AT s» 10K recently observed in experi- 
ments. This is reasonable given the simplified mixture 
model, difference in geometry and many approximations 
of unknown quantities, e.g. a, A7, and S. Also, in ex- 
periments the asymmetric binodal curve of water-2,6- 
lutidine mixtures is much "flatter" near T c compared to 
the binodal in the symmetric Flory-Huggins model we 
use. 

Better agreement with experiments is achieved when 
the values of n£,b are compared. Attraction is observed 
experimentally when is in the range K^b ~ 0.85 — 1 
for cf>Q = 4> c [7] . In preliminary experimental data [12] we 
find K^b s» 0.65 — 0.8 for <j) < (p c and ~ 1.1 — 1.8 for 
4>o > 4>c- These values are in the range of those in Fig. 
[2] and Fig. [3j The difference in temperature range be- 
tween the experiments and our theory can be explained 
by the different scaling of £j = £0 l T l • The experimen- 
tal value of £0 and the critical exponent v = 0.61 [7 
are quite different from the mean field values of £0 and 
v = 0.5. Thus, the adsorbed liquid layers in experiments 
are thicker, leading to attraction at temperatures higher 
above T c . Nevertheless, our results qualitatively show 
that including ion solvation in the theory accounts for 
colloid-surface attraction far from T c in salt containing 
mixtures. 



C. Interaction between hydrophilic and 
hydrophobic colloids 

Another physically relevant scenario is that of hy- 
drophilic and hydrophobic colloids, i.e antisymmetric 
short-range chemical boundary conditions. For this case 
the adsorption force in the absence of salt is repulsive [321 - 
[34] . However, experimentally it was observed that in the 
presence of salt the interaction potential becomes attrac- 
tive when the temperature is decreased toward T c , but 
then repulsive again as T c is further approached [7]. In 
order to explain this phenomenon, Pousaneh and Ciach 
[TJ] assumed ions that are insoluble in the cosolvent and 
suggested that the attraction is due to the hydration of 
ions and is of entropic origin. Another explanation was 
given by Bier et. al. [T3], who attributed the attraction 
to the difference in preferential solvation of ions Au d . In 
their theory the attraction is linear in Au d . 

Fig. [5] (a) shows the non-monotonous behavior of 
U(D) for antisymmetric surface fields within our theo- 
retical framework. Here we assumed the ions have dif- 
ferent solubilities, Au d = —4, and the surface charge of 
the hydrophobic surface is much smaller than that of the 
hydrophilic surface, as is usually the case. For these pa- 
rameters linear theory, relying on the asymmetry in the 
solvation energy of ions, reproduces surprisingly well the 
non- monotonous behavior [13j . In the nonlinear regime 
Au d 7^ is a requisite for this trend but the interaction 
depends on the nominal values of Au* as is shown below. 
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FIG. 5. Inter colloid potentials U(D) for hydrophilic and hy- 
drophobic colloids (antisymmetric boundary conditions). For 
the surface on the right we used Ajr = O.lT/a 2 and ctr = 
— a sat . For the surface on the left we used A^l — —OAT /a 2 
and (jl = —0.01a 3a t <C o aa t- (a) The interaction potential 
U (D) at different temperatures r showing that U becomes 
attractive when r decreases, but repulsive close to T c . Here 
we took for the ions Au + = 4 and Au~ — 8. (b) U{D) at 
r = 0.048 and different values of Au^ . The interaction is 
purely repulsive for Am = Am + — Au~ = (dash-dot curve) 
and weakly attractive for Au d = 2 (dashed curve). The at- 
traction is much stronger in the solid curves, all having differ- 



ent values of Au* but the same difference Au d 



-4. Among 



these curves, the attraction is strongest for the antagonistic 
salt (Au~ = -Au + = 2). 



The onset of the repulsive force is at a surface distance 
D 4£t>, occurring when the two adsorbed solvent layers 
begin to overlap |35j . At larger surface separations, the 
depletion of the more hydrophilic anions close to the hy- 
drophobic surface gives rise also to an enhanced and en- 
ergetically favorable hydrophilic solvent depletion. The 
net result is an attractive interaction at D > 4£b- As the 
critical temperature is approached and ^ increases the 
repulsive adsorption force dominates this relatively small 
attraction. 

Fig. [5] (b) demonstrates the influence of the value of 
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Au + , Au~ , and their difference Au d on the interaction 
profile U(D). In the solid blue curve weplot U(D) for 
Au d = —4 where Au + — 4 as in Fig. pi (a). Reduc- 
tion of Au d from —4 to —2 (dashed curve) diminishes 
the strength of interaction by an order of magnitude, a 
strong nonlinear response. For Au d — (dash-dot curve) 
the interaction is repulsive. Fig. [5] (b) shows that the 
nominal values of Au + and Au~ are important, not just 
their difference. The solid curves in this figure all have 
Au d = —4 but different values of Au^. For a more hy- 
drophilic salt (purple curve, Au + — 6) the attraction is 
weaker while for an antagonistic salt where the cations 
are hydrophobic and the anions are hydrophilic (yellow 
curve, Au + = —2) the attraction is much stronger. The 
reason for this is that hydrophobic cations reduce the wa- 
ter adsorption on the highly charged hydrophilic surface, 
thus reducing the repulsive adsorption force and amplify- 
ing the asymmetric solvation effect. The antisymmetric 
case is an example of the delicate and complex interplay 
between solvent adsorption and electrostatic interactions 
in confined salty mixtures. 



D. Ion specific effects 

In this section we discuss some consequences of the 
specific nature of the ion solvation. The ion solvation 
energy can vary greatly depending on the chemical na- 
ture of the ion and solvent, its value being typically in 
the range Au ~ 1 — 10T. Thus, the influence of the ion- 
solvent coupling on the interaction potential is highly ion- 
specific. In Fig. [6] (a) we plot U (D) for hydrophilic ions 
and surfaces, Au^ = 2, Au~ = 4 and Ajr.l = 0.1T/a 2 . 
When the surfaces are both positively charged (dash-dot 
curve), the anions are in excess between the plates and 
U(D) is attractive since Au~ = 4 is large enough and the 
solvation-related force overcomes electrostatic repulsion. 
However, this is not the case for negatively charged sur- 
faces (dashed curve) where U(D) is repulsive. In this case 
the cations are in excess between the plates and Au + = 2 
is not large enough to overcome the repulsion. Thus, the 
difference in the preferential solvation of cations and an- 
ions can produce a selective interaction with respect to 
the sign of the surface charge. 

An even more intriguing scenario is that of oppositely 
charged surfaces, where for hydrophilic surfaces the inter- 
action is expected to be purely attractive since both ad- 
sorption and electrostatic forces are attractive. Nonethe- 
less, if the mixture contains a strongly antagonistic salt 
(Au + Au~~ < 0, \Au d \ > 1), U(D) has a repulsive re- 
gion, as is seen in the solid curve in Fig. [6] (a) for which 
Au~ = — Au + = 6. In this curve, since <tr = —<jl the 
adsorption of both anions and cations is significant and 
similar in magnitude. Hence, since Au~ = — Au + the 
contribution of the solvation energy to the interaction is 
small. In addition, at separations D > 4£t,, prior to the 
merging of the adsorbed solvent layers [35] , there is sig- 
nificant solvent depletion near a positively charged wall 
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FIG. 6. (a) The effect of the sign of the colloids' charge on the 
inter-colloid potential U(D). The interaction is attractive for 
two positively charged surfaces and is repulsive for two neg- 
atively charged surfaces; compare the dash-dot and dashed 
curves. We used Au + = 2, Ait" = 4, Ay R , L = O.lT/a 2 
and r = 0.008. In the solid curve the surfaces are both hy- 
drophilic (Ajn^L = O.lT/a 2 ) but oppositely charged. For an 
antagonistic salt with Aw - = —Au + = 6 there is a repulsive 
regime at an intermediate range, (b) The effect of the hy- 
drophilicty or hydrophobicity of the colloid's surface. In the 
absence of preferential solvation (Au = 0) two hydrophobic 
(and hydrophilic, not shown) surfaces weakly attract (dash- 
dot curve). For a hydrophilic salt (Ai^ = 4), hydrophobic 
surfaces repel (solid curve) whereas hydrophilic surfaces at- 
tract (dashed curve). We used r = 0.003 and ctl,r = —o~ sa t- 



due to hydrophobic anions. This depletion reduces the 
adsorption force and the overall result is a repulsive inter- 
action. At small separations the adsorption layers merge, 
r increases, and combined with the electrostatic attrac- 
tion the potential is strongly attractive as expected. 

Similar to the effect of the sign of the surface charge, 
the ion-solvent coupling renders the interaction specific 
also relative to the sign of the surface field. In the ab- 
sence of salt, if both surfaces are hydrophobic or hy- 
drophilic they attract [51Tl34j . As we saw before, for hy- 
drophilic surfaces the preferential solvation of hydrophilic 



ions greatly enhances the interaction far from T c , see 
the dashed curve in Fig. [6] (b) for Au* = 4 and 
&7r,L = O.lT/a 2 . By reversing the sign of A 7^ [solid 
curve in Fig. [6] (b)] the interaction becomes repulsive for 
hydrophobic surfaces and hydrophilic ions. The dash-dot 
curve in Fig. [6] (b) shows that for hydrophobic surfaces 
in the absence of preferential solvation (Air* 1 = 0) the 
interaction is weakly attractive. Here, the addition of 
hydrophilic ions renders the interaction repulsive by re- 
ducing the adsorption of the solvent on the hydrophobic 
surface while greatly increasing the adsorption and at- 
traction for hydrophilic surfaces. 

IV. CONCLUSIONS 

In summary, we calculated the interaction potential 
U(D) between two charged colloids in salt-containing bi- 
nary mixtures. The preferential adsorption of one of the 
solvents on the colloid's surface combined with the pref- 
erential solvation of the ions gives rise to a strong attrac- 
tion between the colloids far from the coexistence curve 
of the mixture. For surfaces with symmetric chemical 
affinity, the interaction is governed by the individual ion 
solvation (oc Au*) and to a lesser extent by the differ- 
ence in solvation energies of cations and anions. This 
is true also for the attraction for antisymmetric bound- 
ary conditions albeit here Au d ^ is a requisite for an 
attractive potential. 

The ion and solvent densities near the surface of the 
colloid are highly sensitive to the ion solvation energy. 
Thus, the bulk mixture composition 4>q and salt con- 
centration uq play an important role in determining the 
inter-colloid potential (Fig. [3] and Fig. |4| and this 
suggests the possibility of fine-tuning the potential with 
these readily controllable parameters. In addition, the 
ion specific nature of the solvation energy renders the 
interaction potential sensitive to the sign of the surface 
charge or surface fields, as is shown in Fig. [6j We show 
that the interplay between surface and ion-specific inter- 
actions may lead to non trivial effects whereby (i) two 



repulsive interactions combine to an attraction [Fig. [2] 
(a)], (ii) for a hydrophilic salt, hydrophobic surfaces re- 
pel whereas hydrophilic surfaces attract [Fig. |6](b)], (hi) 
two oppositely charged colloids repel down to distances 
~ lOnm and have a high repulsive barrier of 5-10T [Fig. 
[6] (a)], or (iv) two similarly charged colloids repel when 
they are negatively charged and attract when they are 
both positive [Fig. [6] (a)]. 

Our mean-field theory allows the semi-quantitative in- 
terpretation of recent experiments performed in salty 
mixtures not too close to the critical temperature [7J [T^] . 
We do not look at the influence of the critical fluctuations 
on the interaction between surfaces [6] 13TH34] although it 
has been argued that the addition of salt does not alter 
the universal critical behavior of the solvent [T5] , 

Recently, multilamellar structures were observed in a 
bulk binary mixture upon the addition of an antagonis- 
tic salt [35]. In light of the current work, we believe that 
the effect of antagonistic salts on the interaction between 
surfaces is intriguing and deserves similar attention. In 
addition, for temperatures below the critical temperature 
and/or near the wetting and prewetting transitions, first 
order capillary condensation [10] and bridging transitions 
[11] have been predicted. In light of the current work and 
these recent findings we stress the importance of prefer- 
ential solvation in salty mixtures. Since the Gibbs trans- 
fer energy of ions is usually larger than the thermal en- 
ergy, the solvation-related force can cither induce large 
composition perturbations by itself [10) or amplify them 
significantly if already present, as in this work. 
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